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SUMMARY 


The flow over infinite wedges is investigated theoretically to test 
a conclusion previously reached hy fhe use of Schlichting’s approximate 
method for the caIc\ilation of the laminar boundary layer; namely, that 
in a region of falling pressiore a thick velocity profile can he more 
stable than a thin profile although the velocity at the edge of the 
boundary layer and the pressure gradient are the same for both profiles. 
By the use of the velocity profiles obtained "by Hartree from a numerical 
solution of the boundary-layer equations for wedge flows and by the use 
of Lin’s rapid method for the calculation of the critical Reynolds num- 
ber of a velocity profile, the result is obtained that a thick velocity 
profile on one wedge can be more stable than a thinner profile on a 
wedge of different angle although the velocity outside the boundary 
layer and the pressure gradient are the same for both profiles. This 
result agrees with, and hence confinas, a concltision reached by the -use 
of Schlichting’s approximate method. 

The investigation also leads to the inference that the calculated 
effects of a change in boundary-layer thickness on the stability and 
on the local rotighness Reynolds nu^er should be essentially unchanged 
by replacing the Schlichtlng single-parameter family of velocity pro- 
files by the Hartree e ingle -parameter family of velocity profiles. 


INTRODUCTION 


Experimental investigations have shown that large extents of laminar 
flow can be obtained on airfoils by making the pressure decrease along the 
surface in the direction of flow, by drawing part of the boundary layer 
into the airfoil interior either through a porous airfoil surface or 
through slots cut into the surface, or "by a combination of these methods. 
The experiments, however, make it clear that, to avoid early transition 
to turbulent flow, the siirface must be free of noticeable roughness par- 
ticles and other departures from smoothness. 
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The Indications axe that the dlsturhlng effect of a roughness par- 
ticle can he decreased hy an Increase In the thickness of the boundary 
layer. An Increase In houndaiy-layer thickness, however. Increases the 
houndary-layer Reynolds number and. If the shape of the velocity profile 
Is tinchanged, decreases the ratio of the local critical Reynolds number 
to the local houndary-layer Reynolds number. The assung)tlon Is made that 
the ratio of the local critical Reynolds number (a function of the shape 
of the velocity profile only) to the local houndary-layer Reynolds num- 
ber Is a meas\ire of the stability; the stability decreases as this ratio 
decreases. For a fixed velocity-profile shape, the stability of the 
boundary layer th\is decreases as the thickness Increases. 

The shape of the velocity profile and the thickness of the boundary 
layer at a point on a porous or nonporoxis surface with an arbitrary pres- 
sure distribution can be calculated approximately by the Schlichtlng 
method (ref. l) . This method states that on an inpervloxis surface the 
shape of the velocity profile is determined by the local effective pres- 
s\ire gradient; this gradient is directly proportional to the product of 
the actvial pressure gradient and the square of the boundary-layer thick- 
ness. It is apparent, therefore, that in a region of falling pressure 
the Schlichtlng method predicts that an increase in boundary-layer thick- 
ness increases the effective pressure gradient and thus results in a 
more convex velocity profile. Because the increase in convexity is known 
to Imply an increase in the critical boundary-layer Reynolds number, it 
wo-uld appear that an Increase in boundary-layer thickness could increase 
the local critical Reynolds niomber more than the locsil boundary-layer 
Reynolds mmiber. 

The computations of reference 2, which were made to investigate 
this possibility, led to the resvilt that where the pressure decreases 
along an impermeable surface an increase in boundary-layer thickness 
at a point can, because of the resulting increase in effective pressiure 
gradient, change the velocity-profile shape enough to Increase the ratio 
of the local critical Reynolds number to the local boxmdary-layer Reynolds 
number, emd hence Increase the stability. 

Because the results of reference 2 follow from the Schlichtlng method, 
an approximate method, it is desirable to examine the available exact 
solutions of the boundary-layer equations to determine whether solutions 
exist that can be used to test the results of reference 2. The exact 
solutions of the boundary-layer equations, however, are few, and the only 
solutions that were found suitable it)r tise as a test were the solutions 
for the flow over infinite wedges. 

The method by which the wedge flows are used to test the results of 
reference 2 makes use of the fact that the Schlichtlng method states that 
the shape of a velocity profile is completely determined by the local 
value of the nondimens ional pressure gradient. The previous history of 
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the boundary layer enters only Indirectly. Therefore, the result of 
reference 2 can also he stated in a different manner. Thus, consider 
two inconpresslhle flows with the same kinematic viscosity and let one 
flow contain a point such that the velocity outside the boundary layer 
and the pressure gradient along the surface are equal, respectively, to 
the velocity outside the boundary layer and to the pressure gradient at 
a point in another flow. Then, if the boxmdary-layer thickness is not 
the same at the two points, the thicker boundary- layer velocity pTOfile 
can be the more stable one. 

In the present investigation a confirmation of this conclusion is 
sought by examining the exact solutions of the boundary-layer equations 
for wedge flows, the special set of flows in which the velocity outside 
the boundary layer varies as the distance from the stagnation point 
raised to a power. Numerical solutions of the boundary-layer eq.uations 
for these flows have been given by Hartree (ref. 5 ) • 

Whether the resiilts of reference 2 are always valid is not 
tested in the present work because the wedge flows satisfy a basic 
asstmrptlon of the Schlichting method; namely, that the velocity profile 
is a single-valued fimction of the local effective pressure gradient. 
Therefore, if in ceiisain cases the rate of change of the velocity pro- 
file along the surface would make inaccurate the assun^stion that the 
velocity-profile shape depends only on the local effective pressxnre 
gradient, significant errors might be introduced into the results of 
reference 2. These errors cannot be uncovered by the use of the wedge 
flows. The wedge flows consequently provide only a partial test of 
the results of reference 2. 

In the present Investigation information is also obtained concerning 
the effect on the conclusions of reference 2 of replacing the Schlichting 
single-parameter family of velocity profiles by the Hartree single- 
parameter family of velocity profiles. This portion of the Investigation 
is desirable becaizse the critical Reynolds number is sensitive to the 
shape of the velocity profile. An accurate mnnerlcal calculation, based 
on Lin’s rapid method, is made of the critical Reynolds numbers of all 
the Haxtree velocity profiles. The present investigation is restricted 
to the case of no flow through the surface. 


SYMBOLS 


a 


_g_ 

— 2-p 
ac ^ 


U. 


constant, equal n\imerically to velocity at x = 1 


a 



if 
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coefficients in expression for g 

- 

c 

reference length 


F 

function of Y 


g = J- = < 

+ C^K + CgK^ 


h 

height of roughness particle 


k = l!2 

_ 1 dU 


V dx 

Be U2 dx 


II 

1 

«U'= 



K 

Schlichtlng velocity-profile shape parameter 


u 

velocity in direction of x 


U 

velocity at outer edge of houndary layer 

* 

u = -i 


m 

Uo 



Uo 

reference velocity (velocity at distance c from apex of 
reference wedge) 



velocity of disturhance in boundary layer 



velocity at y = h with roughness particle absent 



reference Beynolds number, UqC^v 


Bg =1 




critical Beynolds number, value of Bq at which a small 



distxirbance is neither damped nor amplified 

- 
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Hh 



c 


y 

y 


roughness Reynolds number. 



velocity through surfa,ce, positive outward 

distance eilong surface, measured from stagnation point 


distance normal to surface, positive outward 


nond lm ens lonal distance from surface. 


y 

)|2 - p 



P 

Po 

8 * 

Si 


parameter 

value of parameter p for reference wedge 


displacement thickness. 


f 

'^0 V 


|(3y 


a measure of hoimdaxy-layer thickness 


e 

V 

a 


momentum thickness of boundary layer. 



function of Y 

stream function 

kinematic viscosity 

angle between walls of wedge, radians 


Subscript: 

w at surface 

Primes denote differentiation with respect to Y; barred quantities 
axe dimensional. 
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ANALYSIS 

Con5)arison of Wedge Flows 


The first portion of ttie present investigation is restricted to 
those inconpressfble laminar flows for which the velocity outside the 
hound ary layer is related to the distance from the stagnation point hy 
the equation 


U 



= ax 


( 1 ) 


which is obtained from reference h, pages I59 to l42, hut put in the 
notation of the present paper. There is no flow throu^ the surface 
in the present analysis. For p ^ 0 the flow is that over the inner 

or outer face of an infinite wedge with an angle 2o! between the walls 
(fig. 1 (a)). The relation between cc and the parameter p is 


a = 

2 

This relation also occurs, in a different notation, on page 5 of refer- 
ence 5. Only values of p less than 2 have physical significance. For 
P < 0 (fig. l(b)) the flow begins at a value of x greater thEin zero 
with a bo\mdary-layer-proflle shape that depends on the parameter p. 

In the wedge flow_defined by equation (l) the velocity U and the 

velocity derivative — are single-valued ftmctions of x. Consequently, 

dx 

two points with equal t) and ^ cannot be found on the same wedge. 

dx 

If, however, two wedges with different angles are considered, and. if the 
velocity at a point on one wedge is properly related to the velocity at 
a point on the other wedge, then two points can be found, one on each 

wedge, with the same U and 

dx 

In order to find the relation between the wedge angles and the two 

points, one on each wedge, that have equal U and — , let 

dx 




(2) 
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and 



(5) 


where U- 


'X2 



are the nondlmens lonal velocity and velocity 


are 


gradient, respectively, in one flow at x = Xq, and and 

the same quantities in the other flow at x = The reference veloc- 

ity Uq and the reference length c are the same in hoth flows. The 
reference velocity Uq can he chosen as the velocity at a distance c 
from the apex of a fixed reference wedge of Incltided angle The 

quantity a (eq. (l)) has the value Sq for the reference flow Eind it is 
found from Uq, c, and Pq hy the use of equation (l) ; 



-2-Po 

c 


The quantities Uq, c, and Pq are arbitrary hut when once chosen are 
fixed. 

For p = Pj^, equation (l) becomes 


- _ 2-Pt 
"• 


or 



Pi. 


2-Pt 

L*1 


(h) 
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where 


U- - 

Xj^=C 


ai = — 


Th 


For p = P 2 ^ equation (l) 'becomes 




P2 

2-P2 


(5) 


where 


%=c 

S2 = -n 


P 2 

_ - 2 -P 2 

E2C 


The condition that U at he equal to U at 
expressed hy equating equations (if-) and (5): 


Pi 


P2 


2-61 2 -Bp 

a^Xi PJ- = agXg 


( 6 ) 


The requirement that at x-i he equal to at X2 (eq. (5)) 


dU 


dx 


expressed hy calculating 


dx 


(-) 

' W/x^ 



from equation ( 4 ) and equating it to 


'xi 


calculated from equation (5) • The result is 


P. 


2 - P 2 


g(P 2 -l) 

'2 2-p2 


Pi 


= a.-, 


2 - 






2 (Pi-l) 

2 -Pi 


(t) 
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When eq,uatlon (7) Is divided by equation (6) the resialt is 


^2 _ P2 ^ - Pi 
Xi Pi 2 - p2 


(S) 


where, in order for equation (5) to be satisfied, p^ anri pg must be 


either both positive or both negative. 


Thus, U and — at x = Xo. 
dx 


in a flnw with p = P2 and a = & 2 , are equal, respectively, to U 
dx 


and — at 


X, = 


Pi 2 - p 2 


^ r- 


P 2 2 - Pi 


in a flow with p = Pi and a = ai- The relation between ai and a^ 

is obtained by the use of equation ( 8 ) and either equation ( 6 ) or equa- 
tion (7) . The result is 


ag - aiXi 


2 (Pi-P 2 ) 

(2-Pi) (2-P2) /^ ^ ~ ^2 

\P2 2 - Piy 


P 2 

2 -P 2 


( 9 ) 


Consequently, in the flow over infinite wedges a point on one wedge 

can have U and ^ equal, respectively, to U and ^ at a point on 
dx dx 

a wedge of different angle. The arbitrary quantities ai, Xi, Pi, 

and Pg determine Xg and ag by equations ( 8 ) and (9). Fixing 

the arbitrary qu ant ity ai is eqxiivalent to fixing u/Uq at x = c 

in the flow for which p = Pi* The quantity eig is not arbitrary but 

is determined by equation (9); the quantity ag is the value of 
U^Uq at X = c on the wedge with p = pg. 
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Stability of 


-Boxmdary-Lajrer 


Velocity Profiles on Wedges 


The ass\jnrption concerning the stability of a 1 ami nar-bonnd ary-layer 
Telocity profile is the same as in reference 2, namely, that the ratio 
Rq jHQ at a distance x from the stagnation point is a measure of the 


stability of the boundary layer at that value of x and that the boundary 
layer becomes more stable as the ratio Rq ^Rq Increases. In order to 

congoare the stability of a velocity profile at x = Xg with the sta- 
bility of a profile at x = Xj^, it is necessary to find Rq eind Rq 


at both Xj^ and Xg. 


In order to obtain Rg 
tion is 


the definition for 0 is xisedj the defini- 


6 



The velocity u 


is obtained from the relations 


and 



V = - P 


1 

_l/2 5 Ib 1/2 
a X V F 


y = ^2 - p 


i=P 

_-l/2_2-B_l/2 
a ' X V ' y 


( 10 ) 


( 11 ) 


vhlch 8ire from reference h, pages 159 to 142, but the notation has been 
cheinged to that of the present paper. 
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The result is 


a = Sx^"^' (12) 

where the function F(y) satisfies the differential equation (ref. k-, 
pp. 139 to 142): 


Jim ^ jrptt 


- p[(F')^ •* 3 


0 


(13) 


and the houndary conditions 


Y = 0 


F' = 0, F = 0 


F' 


When p = 0, eq.uation (13) becomes the Blasius equation for the flow 
over a flat plate (ref. 4, p. 135) • When equation (l2) is used with 
eq.xiation (l) in dimensional form. 


_g_ 

,-T — 2-p 

U = ax ^ 


(14) 


the result is 



(15) 


By making use of eq.mtions (ll), (l4), and (l5)^ the expression for 0 
can he written as 


i\ = ^2 - p r F'(l - F’)<ai 

X y V 'Jq 


(16) 
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The expression for Eq then is 


V 


^ ^2 - P r F‘(l - F')dj 

Z J r\ 


or 


Eg = ^ ^2 - p J F*(l - F’)< 3 y 

From equation (l) it follows that 


(17) 


2 

The = 6oc^“P 


and that 

IS. 2 

° ^ ^ 
^ g(P-l) P 

dx p 2-p 

a — t- — X 


Therefore : 

Ux = — 2 — j£. 

2 - p ^ 

dx 

Expression (I7) for Eg can thus he written as 



(18) 
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The integral / F’(l - a function of was calculated from 

table II of reference 3 1^7 ase of Sin5)Son's rule (ref. 6 , pp. 120 
to 122) . All the entries in the table for a specific value of p were 

used to get the integral for that value of p. The function ^ 


dx 

I • w 

is shown in figure 2 and the integral / F'(l - F‘)<H is given in 


TJ2r, 


table I. The integral / (l - F')dl, related to the displacement 


0 

thickness 5* by the eq.uation 


= ^2 - p r” (1 - F')< 3 Y 

X 1 V <J0 


is also given in table I. 

In order to find the expression for the conclusion of refer- 

ence 7 that the critical Reynolds number of a velocity profile depends 
only on the shape of the velocity profile is used. For wedge flows the 
shape of the velocity profile depends only on the value of the param- 
eter p (see eq. (13)). Therefore, the critical Reynolds number Rq^ 

depends only on p: 



The critical Reynolds number is calculated by Lin's rapid method 
(ref. 8) . The critical Reynolds ntmibers for some of the velocity pro- 
files given by Hartree in reference 5 have already been calculated by 
Lin's rapid method and are given in figure 2 of reference 9- Because 
only the first significant figure of the critical Reynolds number can 
be read without estimation from this figure and because critical 
Reynolds numbers were calculated for less than half of the profiles 
given by Bartree in reference the critical Reynolds numbers eire com- 
puted in the present work for all the Bartree profiles. 
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Lin's sinplifled formula for Eq is 



\ii/ 


(19) 


wMch. is taken from reference 8 tut has "been pvrt in the notation of the 
present paper. 

This formula is to he put in a form that contains the VEiriahles F 
eind Y of equation (13) • The use of eq[uations (ll), (l^), (l 5 )^ (16) 

leads to 


/at' 


u — \ pco 

— 2 . = F« / F'(l - F')dY 

\^ll ^0 


(20) 


■ The term ^ eq.xiation (19) is the value of u/u for which 

the relation 


-It 


/a|\ 

u 




5 - 




i 

2 

— — 

y 



1 9 

\ 0/ 

w 

u 


U 



u 


82 b 
u 





= 0.58 


(from ref. 8 ) is satisfied. 

From eq.\mtions (u), (l 4 ), and (I6) the relation 

Y 


0 


j poo 

F'(l - F')dY 

0 


(21) 


(22) 


is obtained. 
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By using eq^uations (20) and (22) the relation 



is obtained. 


The equation for Ucju (eq. (2l)) can now be written 


as 




^w"^F'F"' 

= 0.58 

F’ y(F")5 


(25) 


and the equation for R0^ (eq. (19)) can be bitten as 



where F^' is the vaJLue of F' for which equation (25) is satisfied. 

The values of Fv" given in table II were taken from Hartree 
(table II of ref. where Fw" is called y" ( O) ) . The integral 

/ F'(l-F')<i7 has alreaidy been con 5 )uted (see eq. ( 18 ) and table l) . 

A detailed description of the computation of F^,' is given in the 
appendix. For p ^ -0.l4, numerical differentiation of Hartree 's tab- 
ular values of F' (y' in Hartree ' s notation) was not used to find F" 
and F'” in the equation for F^' (eq. ( 25 )). - Instead, the func- 
tion F' was expanded from Y = 0 in a Taylor's series and the deriva- 
tives of F ' were found by differentiating the series . This procedure 
is believed to result in more accurate values of F^.', apd hence of 

than could have been obtained by numerically differentiating Hartree ' s 
tabular values to find F" and F .. The variation of R 0 ^ with p, 

calculated by the equation for (eq. (24)), is shown in figure 5- 

The values of R 0 ^ are also given in table H. 
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The ratio B0 ^|Eq can now be found by dividing equation (24) by 


equation (18); the result, — ^ 1 

% \ 


U^Rf 


dx 


, is shown as a function of p in 


figinre 4. Figures 2 anfl 4 can be used to determine whether the flows 
defined by equation (l) contain points with equal U and — and with 


the property that the point with the larger value of 0 has the larger 
value of 'RQ^j'RQ'f the value of is the same for all the flows. Flg- 

■ure 2 shows that when two points have equal values of R^, U, and ^ 


but different values of p, the point that has the larger absolute value 
of p has the larger value of Rg. It is recalled that because of 

equation (3) both values of p have the same sign. Figure 4 shows that 
for points with the same values of R^, U, and ^ the value of RQ^^Rg 


increases with p for p > O.O5. Therefore, for p > O.O5 the point 
with the larger value of Rg has the larger value of R0^Rgj thus, of 

two points with equal values of U and — in flows with equal R^, the 

dx 

point with the larger value of 0 can have the more stable velocity pro- 
file, The exact solutions of the boundary-layer equations for a special 
type of pressure distribution thus confirm a result of reference 2. 


Con^jarison of Predictions of a Method Based on the Hartree 
Profiles With Predictions of the Schlichting Method 


The purpose of the second part of the present investigation is to 
determine whether replacing the Schlichting single -parameter family of 
velocity profiles (ref. l) by the Hartree single-parameter family of 
velocity profiles (ref. 3) would leave the conclizsions of reference 2 
essentially unchanged. 


The investigation is ne.de by first deteiminlng whether or not an 
approximate method based on the Hartree velocity profiles predicts that 
an increase in boundary-layer thickness alone at a point can increase 
the stability of the boundary layer at that point. The values of Rq^ 


for the Hartree profiles are then compared with the values of Rq^, for 

^ R0C fu^ 

the corresponding Schlichting profiles, and the curve of - — 


dx 


for 


the Hartree profiles is conpared with that for the Schlichting profiles . 
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Of Interest also is a comparison of the rate of change 
profile shape parameter and of the stability parameter 


of the velocity- 

Iu2r, 


1 

1 dx 


with 


cha n ge in boundary- layer thickness. Finally, the rate of change of the 
local roughness Reynolds number with change in boundary-layer thickness 
is calculated for the Hartree prafiles and con5)ared with that for the 
Schlichting profiles. In order to set up a method based on the Bartree 
profiles, only the boundary-layer momentum equation (ref. 4 , p. I35) and 
equation (18) need be used. Equation (18) states that the local value 

of 9 determines the local value of p. Positive values of p 

U%c 

dx 

occior for positive values of — . The determination of p completely 


determines the local state of the boundary layer. A method based on 
the Hartree profiles is given in references 10 and 11 . 


In order to determine whether or not an approximate method based 
on the Hartree velocity profiles predicts that an increase in boundary- 
layer thickness alone at a point can increase the value of Bq^IKq at 

that point, only figures 2 and It- need be used. In figure 2 is shown the 


connection between p and 


Rc 


M 

dx 


the figure shows that p increases 


as 


Re 




increases when p is positive. Therefore, p increases 


dU 

Hx 


as 0 increases, when p is positive and R^, U, and 

U^R^ 


dx 


are fixed. 


Re 

Figure 4 shows how the stability pareimeter — - 

Re 


for p > 0.05 an increase in p Increases 


R0c 




dxl 


dx 


depends on pj 


Hence, an 


increase in 0 , with R^, U, and ^ fixed, results in an increase 


dx 


in p, and consequently, for p > O.O5, in an Increase in 'Rq^I'Rq. 


On 
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the other hand, the combination of figures 2 and 4 shows that when 
0 < p < 0.05 a small Increase in 0 results in a decrease in R0^^R0. 

This behavior is the same as that found in reference 2. 


In order to con^jare the ciirves for R0^ 


fl-Tirl 


^9c 




dx 


obtained 


in reference 2 as functions of the 
eter K with the CTxrves for R0^ 


Schlichtl ng ve locity-profile param- 

R0C Ju^Rc 


and 


'0 


llil 


obtained in the pres- 


ent investigation as functions of p, both K and p are replaced by 
a parameter called k by Schlichtlng (ref. l) and k in reference 2. 
The parameter k is the value of the nondimens ional pressure gradient 
(ref. 2): 


k = ^ ^ 

Rc xj^ 


(25) 


The relation between p and k can be obtained by combining equa- 
tions (18) and (25); the result is 


k = p 



The relation between K 


and k is (ref. 2) 


(26) 


k = g2(K + 1) (ki = 0) (27) 


Therefore, the value of p and the value of K that occvir together for 
any value of k satisfy the equation obtained from equation (26) and 
equation (27); thus. 


g^(K + 1) 


P 



F’)^ 


= 0 
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where 


g = Cq + C^K + C2K^ 


G?he values of C 


C2 are given in reference 1 . The values 
of p and K that occur together are therefore given hy the equation 


and 


(Cq + + C2k2)^(K + 1) - p J f’( 1 - F')dY 


= 0 


( 28 ) 


For each value of p for which Hartree gives a velocity profile (ref. 5) 
the corresponding value of K was found from equation (28) hy a process 
of successive approximation. The variation of p and K with the non- 
dimensional pressirre gradient k is given in table III and in figure 5. 
Figure 6 gives the variation of R0 with k for the Hartree profiles 


Re 




for 


and for the Schlichting profiles. Figure 7 shows how 

dxl 

both the Hartree and the Schlichting profiles depends on k. Because 
the abscissa k in figures 6 and 7 is a function of p alone, fig- 
ures 6 and 7 can also be toterp reted as a comparison between the exact 

fu% 


values of R9 and 


Rc 


of Rq and =-2. 


Re ^ 


Re 


dU 

dx 


for wedges and the approximate values 


dx 


calculated for these wedges by making use of 


Schlichting 's method (ref. l) . 

In order to find how dp/dRg, the rate of change of profile shape 
parameter with change in boundary-layer thickness, depends on p, equa- 
tion (18) is used. The result for Rr, is 

^Rq 


Re 


8Rn 


p 


p 


I 


F'(l - F')*!! 


^ / F'(l - F«)dY + i 
dp Jo 2 


(29) 
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The quantity 



ccilcrulete — f F’(l 

Jq 


F’)dY has already heeii calculated. In order to 
a table of divided differences was formed 


d 

with, p as the Independent variable. The derivative — / F'(l-F’)ciy 

^ Jq 

was foimd from the table by use of the formula of exan5>le k, page 217 of 

reference 6. The function Eg ^ is shown in figure 8. In the same 

dEg 

figvire is shown Eg - — for the Schlichting profiles, taken from refer- 

dEg 

ence 2. Altho-ugh the ratio of Eg to Eg becomes large as k 




departs from zero, the rate of change of 


Be, 




Be y m 

dx 


with Eg (really 0 ) 


is about the same for the Schlichting and Hartree profiles for moat of 
the range of kj thus 



but 


P 2 

iL. i ^ 

Ec dx 


therefore 


^ = 2k 
^ dEg 


(Ec, U, and — fixed) 
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11s 


= 2k 




Because figure 7 shows 


%c |u%c 

FiM 


to he about the same for both the 


Hartree and the Schlichting velocity profiles over most of the range of k, 

it follows from equation (30) that the rate of change of — 

Re ' ^ 

dx 

with 0 is aljBO about the same for the Hartree a-ni^ Schlichtin g profiles. 
In the region nesir k = 0,015, however, there are differences; the deriva- 

./R0^ |d2r,\ 


1 il/ 

tive Rg — is positive for the Hartree profiles for values 

of k for which it is negative for the Schlichting profiles. 

The effect of an increase in boundary-layer thickness on the rou^- 
ness Reynolds nxmiber Rjj^ is obtained by substituting the expression 


V' F'(l - F'] 


for the term f(K) in equation (16) of reference 2. The eq.uation 
Rq ^Bvi 

for -2. --ii for the Hartree profiles becomes, after a development 

^h °Pe 

parallel to that of reference 2, 


% ^h _ 

Rjj SRq 


.1 + Eq ^ 
^ SRe 


f^oo 1 

F^" / F»(i - 

dp ^0 J poo 

F^" / F'(l - F')dl 
^0 
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Figure 9 shows how — ■:: — varies with the nondimens ional presstire gra- 

% SBe 

dient k for hoth the Hartree and the Schlichting velocity profiles; 
the values for the Schlichting profiles are taken from reference 2. The 
use of either the Hartree or the Schlichting velocity profiles leads to 
the conclusion that an increase in houndary-layer thickness decreases the 
roughness Reynolds number. 


The information contained in figures 6 to 9 indicates that a method 
based on the Hartree profiles should predict essentially the same effects 
of an increase in houndary-layer thickness as the calculations of refer- 
ence 2., which made use of the Schlichting method. For example, a calcu- 
lation made hy the Schlichting method predicts that if the value of 0 
S’i 35 percent of the chord from the leading edge of the NACA 64A010 air- 
foil at a Reynolds number of lo'^ is multiplied hy 3*^5^ ilie value of Rq 


is increased from 1,221 to 4,208 and the value of Rq^Rq Is increased 

from 0.239 to 1.0 (see p. 23 of ref. 2). A calculation made hy the use 
of the information for the Hartree profiles in figure 7 predicts that 
this value of 0 shoiild he multiplied hy 3*27 to increase Rq J'Rq 


from 0.210 to 1.0; the value of Rq is Increased from 1,221 to 4,000. 


1 

U%c 

The value of 1 

c 

\ - 


\ 

dx 


at this position of the airfoil for 


Rc = loT 


is 


equal to l4,260. 


DISCUSSION 


The numerlcaJ. solutions of the hoimdary- layer eqmtions for the 
flow over infinite wedges are used to show that a thick velocity profile 
on one wedge can he more stable than a thin profile on another wedge, 
although the velocity outside the boundary layer and the pressure gra- 
dient along the surface are the same for hoth profiles. A result of 
the calculations of reference 2, which were based on the approximate 
Schlichting method, is thus confirmed for a special type of flow. 

It is noted, however, that the wedge flows defined hy equation (l) 
satisfy a basic assumption of the Schlichting method, namely, that the 
velocity-profile shape at ary point on the surface, for no flow throTogh 
the surface, is determined hy the value of k at that point. Equa- 
tion (26) shows that for wedge flows the shape of the velocity profile 



KACA TN 2976 


25 


Is indeed a function of k only. For these flows k is independent 
of X and depends only on the wedge angle; for a fixed wedge angle there 
is thuB no change in k or in velocity profile along x. On the other 
hand, in regions in which k varies rapidly along the surface and in 
which, consequently, the basic assumption of the Schlichting method that 
the velocity-profile shape depends only on k is not satisfied, the 
predictions of the Schlichting method will prohahly he less precise than 
for the wedge flows. For exanple, a discontinuity in the distribution 
of k cannot result in a discontinuous change in the shape of the 
velocity profile. The profile shape predicted by the Schlichting method 
can, however, be ejpected, in a region of falling pressure, greuiually to 
approach a shape that is in the neighborhood of the correct one as the 
distance downstream from the discontinuity in k is increased. 

The result of the second portion of the analysis is that a method 
for the calculation of the houndary layer based on the Hartree profiles 
should indicate essentially the same effects of an increase in boundary- 
layer thickness as the calculations of reference 2 , which were based on 
the Schlichting method. This result makes the concltisions of reference 2 
appear to be less dependent on the particular single-parameter family of 
velocity profiles chosen by Schlichting; one of the weaknesses of the 
einalysis of reference 2 is that it is based on a particular single- 
parameter ffimily of velocity profiles and that Eg is sensitive to 

the shape of a velocity profile. 


CONCLUSIONS 


The result is obtained that in flows over Infinite wedges cases 
occur for which a thick boundaiy-layer velocity profile on one wedge is 
more stable than a thinner profile on a wedge of different angle, although 
the velocity outside the boiindary layer and the pressure gradient are the 
same for both profiles . This result confirms a conclusion previously 
reached by the use of Schlichting ’s approximate analysis, namely, that 
in a region of falling press\ire the stability of a thick velocity profile 
can be greater than that of a thin profile when the velocity at the edge 
of the boundary layer and the pressure gradient are the same for both 
profiles . 

The investigation also leads to the inference that the calculated 
effects of a clmnge in boundary-layer thickness on the stability and on 
the local roughness Reynolds number should be essentially xmchanged by 



ItACA TK 2976 


2k 


r^lacing the Schlichting single-parameter family of velocity profiles 
hy the Hartree single-parameter family of velocity profiles. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
L a n gley Field, Va., May 11, 1953. 
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APPENDIX 

CALCUIATION OF Pg’ AND 

In order to confute Fg* it is necessary first to confute F" 

and F'*' over a range of Y large enough for equation (25) to he 
satisfied at one point in this range. The congjutation of F" and F'" 
was made hy first expanding the function F' in a Taylor's series from 

Y = Oj the series, to Y^5^ is 



(52) 


where F ° means the nth derivative of F at Y = 0. The derivative 
of the series for F' results in the series for F”: 


Ik p iun-2 

F" -Y "" 

^ (n - 2) I 


(55) 


and the derivative of the series for F" results in the series for F' 


^ - 5)- 


(5^) 


The coefficients F^^ sure obtained hy successively differentiating equa 

tion (15) and using the houndaxy conditions at Y = 0. Thus, directly 
from equation (15) and the conditions 


Y = 0 


F = 0 


F' = 0 
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there resrilts: 


= -P 


The first derivative of equation (15) is 


pIV ^ j-pm ^ f'F"( 1 - 2p) = 0 
The "boundary conditions at Y = 0 then lead to 



= 0 


YTV 

This process was repeated over over until was obtained. The 

coefficients axe listed in table TV. The quantity depends only 

on p and is listed in table Hj these values of F^" were taken from 
reference 5 (where the symbol y"(0) was used instead of F^”) . 

The series (52), (35) ^ (5^) were used to calculate F', F", 

and F'” by xising the value of F^.” given by Hartree (ref. 5) for 

each value of p. For p ^ 0, the series (32), (55)^ and {j>k) were 
assumed to be valid representations of F', F", and F"* from Y = 0 

to the largest value of Y at which the magnitude of the last term in 
the series for F"' was less than O.OOOO5. For p = -0.10 the cri- 
terion was relaxed to allow O.OOOO8, and for p = -0.l4 the value O.OOOI5 
was allowed. For values of Y less than the largest values of Y at 
which these requirements were satisfied, the values of F' calculated 
by the series (32) did not show a consistent departure from Hartree *s 
tabulated values (ref. 3) • 

In order to calciilate Y^, the value of Y at which equation (23) 
is satisfied, equations (33) (3^) and the values of F^" and F' 

given by Hartree (ref. 3) were used to calculate the left-hand side of 
equation (23), called 0, for a range of Y large eno^lgh to include 
values of 0 greater than 0.58. For values of -0.l4 < p< 2. it, the 
con^jutation of Y^ and F^,' then consisted of an iteration procedure 

in which Y^, was first estinated by applying a divided-difference 



KACA TN 2976 


27 


method (ref. 6 , ch. 7) to the calculated values of 0. A new value 
of 0 for this value of Yq was con^juted hy calculating F', F", 

and F*" for this value of Y from equations (32), (33), and (34). 

If 0 was not equal, to O.58 the procedure was repeated and tuse was 
made of the new values of 0 and Y in the divided-difference pro- 
cedure. This iteration procedure resulted in values of Y^ for 

which 0 was equal to O.58OO ± 0.0001 for all values of p for which 
-0.l4 -g p g 2.4. The last iteration for Y^ also gave the final value 
of Fc'. 

For p = -O.I6, -O.I8, and -O.19 the Taylor's series is not valid 
to Yj.. In these cases the values of F' in Hartree's table (ref. 3) 

were used to form difference tables. In most cases third differences 
were used. These difference tables were used with Newton's anil Gauss' 
formulas (ref. 6. pp. I60 to I7O) to con5>ute the values of F" and F'". 
The values of F" emd F'" con5)uted by the Newton anil Gauss formulas ■ 
were smoothed by one application of the 7-point formula on page 278 
of reference 6. These smoothed values of F" anil F'" were used with 
Hartree's tabular values for F' (ref. 3) to compute 0. The value 
of Y for 0 = 0.58 was then found by inverse inteip)olation, for which 
a divided-difference formula was used (ref. 6) . For this value of Y, 
called Yg, the value of F', called F^', was foimd by interpolation 

from the difference tables j third differences were used. 

Once Fj, ' was calculated^ the critical Reynolds number was com- 
puted from the equation for Rq^ (eq. (24)). The variation of Rg 

with p is shown in figure 3. The values of Rq^ are also given in 

table II^ together with Hartree's values (ref. 3) of F^". The values 

of Rq^ in table II are believed to be acc\orate to about 6 in 10,000 

for -0.l4 ^ p ^ 2.4. This estimate of the accuracy is based on an 
estimate of the effect on the calculated value of Rq^ of a small dif- 
ference in 0 from O.58OO. For p g -O.16 the values of Rq^ are 

probably accurate to the number of significant figures given in the 
table. These estimates of the acciaracy concern only the numerical pro- 
cedure and do not include the effects on Rg^ of any inexactness in 

Hartree's values of F^" or of any inexactness in Lin's formula (eq. (19)) . 

In the course of the conq)utations it was noted that the value O.8860 
at X = 1.6 (Hartree's notation, ref. 3)' and p = 0.5 did not seem to 
be consistent with neighboring values. An interpolation which made use 
of the other values in the vicinity gave O.876O instead of O.886O. The 
value 0.8760 was used in the computations of the present analysis. 
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TABLE I 


THE INTEGRALS 



- FOiH AND 



F’)< 3 Y FOR 


VARIOUS VALUES OF p 


p 

pco 

/ F’(l - F')< 3 Y 

Poo 

/ (1 - F*)dY 

^0 

2.h 

0.211I-9 

0.4607 

2.0 

.2308 

.4974 

1.6 

.2502 

.5439 

1.2 

.2760 

.6068 

1.0 

.2925 

.6479 

.8 

.3119 

.6987 

.6 

.3358 

.7638 

• 5 

.3502 

.8044 

A 

.3667 

.8526 

.3 

.3858 

.9109 

.2 

.11-080 

.9858 

.1 


1.0803 

0 

.4695 

1.2167 

-.10 

.5155 

1.4436 

-.111- 

.5388 

1.5967 

-.16 

.5523 

1.7076 

-.18 

.5676 

1.8712 

-.19 

.5768 

2.0070 

-.1988 

.5852 

2.5587 
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TABLE H 


THE QUANTITIES F^" AND E 0 ^ FOR VARIOUS VALUES OF p 


p 

V 

(1) 


2.4 

1.857 

6,590 

2.0 

1.687 

6,209 

1.6 

1.521 

5,728 

1.2 

1.556 

5,002 

1.0 

1.2526 

4,529 

.8 

1.120 

5,928 

.6 

.9960 

5,132 

•5 

.9277 

2,662 

.4 

.8542 

2,154 

.3 

.7748 

1,551 

.2 

.6869 

965.0 

.1 

.5870 

470.4 

0 

. 46 ^ 

179.5 

-.10 

.5191 

62.5 

-.14 

.2595 

58.4 

-.16 

.1905 

50 

-.18 

.1285 

17 

-.19 

.086 

7 

-.1988 

0 



^From Hartree's table (ref. 5) 
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TABLE III 

THE VELOCITI-HROFILE PARAMETERS p AND K FOR VARIOUS VALUES 
OF THE PRESSURE-GRADIENT PARAMETER k 


k 


K 

0.1108 

2.4 

- 0.4818 

.1067 

2.0 

-.4986 

.1002 

1.6 

-.5237 

.09142 

1.2 

-.5592 

.08542 

1.0 

-.5839 

.07781 

.8 

-.6159 

.06764 

.6 

-.6597 

.06131 

.5 

-.6877 

.05379 

.4 

-.7217 

.04465 

.3 

-.7644 

.03329 

.2 

-.8196 

.01896 

.1 

-.8934 

0 

0 

-1 

-.02658 

-.10 

-1.176 

-.04064 

-.14 

-1.289 

-.04881 

-.16 

-1.366 

-.05800 

-.18 

-1.470 

-.06321 

-.19 

-1.543 

-.06810 

-.1988 

-1.629 
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TABLE IV 

COEEFICIErraS m tee TAYLQR’B SERIEB for F‘, F”, AKD F'” 



f/ = -(Fi")2(i - 2p) 

F^VI = 2P3^"p(2 - 33) 

Fv^ = -2^(2 - 3p) 

Fy^^ = (Fi") 3(1 - 2p)(ll - lOp) 

F^^ = -2(Fi'’)^p(lt-5 - U3P + 66p2) 

Fv^ = l6F3^"p2(2 - 3P)(8 - 7P) 

F^^ = -(Fi")^( 1 - 2p) [£29 - l6p)(U - lOp) + lif(l - 2p)(4 - 5py] - 
l6p5(8 - 7p)(2 - 3p) 


F^^ = (F3^")5pg(57 - l8p)(l^5 - 115P + 66p2) + (1 - 2p)(U - 10p)(95 - 72p^ 
= -2(F3^")2p2|i6(23 - 10p)(2 - 5P)(8 - 7p) + 

10(13 - 9p)(l<-5 - U-5P + 66p2) + 6(27 - 28p)(l - 2p)(2 - 3P)] 


= (Fi") 5-[2(28 - Hp)(l - 2p)[£29 - l6p)(ll - lOp) + l4(l - 2p)(l). - 5pT| + 
3(139 - 136p)(l - 2p)2(u - 10p)J- + 8F^"p5 [lt(28 - Up) (8 - 7p)(2 - 3p) + 
^(8 - 5P)(2 - 3P)(8 - 7p) + 5(27 - 28p)(2 - 3p)f| 
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Figure 2 .- Variati ons of the houndary-layer Reynolds number 

/ fu^ 

— ^rilth the velocity-profile shape parameter p. 


parameter Eg, 


dUI 
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Figure 5 -- Variation of the velocity-profile shape parameters p and K 
with the nondimensionsLl pressure gradient k. 
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Figure 6 .- Variation of the critical boxzndary- layer Reynolds nvunber Ro 

®c 

irf.th the nondimensional pressure-gradient parameter k for the 
Hartree and the Schlichting velocity profiles . k = 

Rc Ij2 dx 






Hartrae TaXooit; profilaa 
Sohliohtlng Telooitj profllaa 
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Figure 8 .- Variation with the nondimens ional pressure gradient k 

of Rq the parameter for rate of change of p with boundary 

0 8R9" 

layer thickuess, and of Rq , the par a meter for rate of c han ge 

8Rq 

of K with houndary-layer thickness. 
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Hartree velocity profiles 
Schlicbting velocity profiles 


Figure 9«- Variation with the nondimensional pressure gradient k 
Rq SRv, 

of —2. — the parameter for rate of change of roughness Reynolds 

Rh SR0 

number Rj^ with hoimdary-layer thickness, for the Hartree and the 
Schlichting velocity profiles , 


NACA-Lasglev > &>6>>53 • 1000 






